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^ . Abstract: 

The recently observed transient absorption of the lowest excited state of the special pair P* at 
2710c?7z~^ [1] is assigned as a singlet which arises from the coupling of the two lowest triplets from the 
r~| ■ two dimer halves. INDO calculations are used to predict its intensity. The analogy of the coupling 
O \ mechanism to the trip-doublet spectrum from P^ is shown and the influence of the double triplet on 
Y\ ' the Stark effect of P* is investigated. 

>~>. Introduction 

Qh| In recent femtosecond infrared spectroscopy of the reaction center excited states of Rb.sphaeroides up 

to lOOOOcm"^ an interesting spectrum has been observed. A relatively sharp band at 2710c7n~^ is 

followed by a very broad around 5300c7n^^ [1 ]. These bands have been interpreted as internal charge 

transfer states, where the extra broadening of the upper is seen as a result of mixing to higher excited 

states of the dimer. This interpretation can not explain that the higher state has more intensity, 

^^ . also the absence of the lower in the ground state spectrum and in the Stark effect spectrum would 

S^ \ be in conflict with this model, since the CT state should mix appreciably with P* . From our INDO 

f^~. ' calculations which are based on the X-ray structure of rps.viridis we have predicted earlier [ 2 ] the 

^> ■ presence of a low lying excited state denoted P** which is a singlet component of a double excited state 

tZJ \ with a triplet on each of the two dimer halves. We obtained this state too low in energy as we know 

• ^H \ from the predicted energy positions of the lowest triplet. However, if we correct for this shortcoming of 

^.^' the INDO calculation, we predict this state to be at the position of the newly observed band. For the 

,i^ . symmetric dimer it is an even state and therefore does not couple with P* . There exists, however, an 

'^' interesting one-electron coupling mechanism to the internal charge transfer states. This coupling is of 

the same nature as the one which provides intensity to the trip-doublet observed in the P~^ spectrum. 

\^ . In this paper we want to show that the P* — > P** transition obtains its intensity from the CT states. 

H \ Moreover the presence of the P** state helps to understand the unusual Stark effect spectrum of P* 

and provides insight into the internal asymmetry of the dimer. We start with the analysis of the 

spectrum of P"*" , then we present an extended exciton model for the dimer which includes the CT and 

the double-triplet state. Analytical expressions are derived for the infrared transitions and the Stark 

effect. The relevant couplings are discussed in the framework of INDO calculations. 

Simulation of the P+ spectrum 

The spectrum of P^ [1,3] shows a very low energy transition at ~ 2500cm^^ with SODebye^ intensity 
and another transition of similar strength at ~ 7500c?ti^^ with an additional shoulder at higher 
transition energy. INDO calculations by us [2] as well as model studies by Breton et al [3] and by 



Parson et al [4 ] related the lowest transition to the hole transfer. The higher transition shows more 
complex configuration interaction [2]. Recent investigations by Reimers et al [5 ] attributed also the 
low energy transition to the hole transfer between the two dimer halves and the higher band to a 
trip-doublet state. A quantitative understanding is still missing since the internal asymmetries of the 
dimer wich are difficult to assess, influence the relative intensities in a sensitive way. 

Here we want to discuss the experimental spectrum in the framework of a simplified four state 
model which incorporates the two cationic states L"*" and M"*" and the two lowest excitations for each 
of them which are the triplet-coupled doublets L~^M^ and M^Lr". A more detailed analysis including 
higher excited configurations explicitly will be given in a separate paper . 
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Within this basis the positive charge is fully localized on P^ for the first and third state and on 
Pm for the other two. The delocalization of the charge is the result of the intermolecular resonance 
interaction Ui ,U2 competing with the small energy differences 2di^2 of the symmetry related states. 
Such asymmetries arise from asymmetric interaction with the surrounding, e.g. hydrogen bonding to 
protein residues or from asymmetries of the two chromophores. We investigated especially the influence 
of the orientation of the two acetyl groups. A structure model where the oxygen and the methyl group 
of Pm are exchanged is also consistent with the X-ray structure and leads to a large splitting of the 
intradimer charge transfer states which is needed to explain the strong dipole change of P* seen in 
the Stark spectrum. The two couples of states denoted by the indices 1 and 2 respectively interact 
via the matrixelements U12 and U2i- A detailed analysis of the INDO results shows that the coupling 
Ui is mainly given by the negative resonance interaction —(3lm involving the molecular orbitals of 
two different charged states P^ and P^. For our analysis we want to interpret couplings for P^ and 
P* within one set of molecular orbitals. We use localized SCF orbitals of the neutral dimer which 
are obtained switching off all resonance interactions between the two dimer halves during the SCF 
step. Electron correlation effects are important and the expression for Ui becomes more complicated 
as configuration mixing takes place. We find that Ui is underestimated by about a factor of 2 if the 
resonance integral I3lm of the SCF orbitals of P is used. For the coupling of the two trip-doublet 
states the reorganization effect is of little importance. The couplings of the lowest cationic states to 
the trip-doublets involve the transfer of an electron from the HOMO on one side to the LUMO on the 
other side. Here the reorganization effect leads to a significant increase of U12 and U2i- These rather 
strong couplings induce the intensity of the trip-doublets. The results are summarized in table 1 
For a simplified analytical treatment we first diagonalize the two pairs of states with their asymmetries 
which we want to determine empirically. The couplings C/12 and U21 are treated as a perturbation. 
We use the linear combinations of states 
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Specifically one obtains (pj = O(ib^) if the lower state of a pair is localized on Pum) &iid 'Pj = if if 
both states of the pair are fully delocalized. Within this model we can determine the main parameters 
empirically and compare them with our calculated ones. From the experimentally observed charge 
distribution of « 65% on Pl we get sin^((/>i) f» 0.35. 

The intensity of the lowest transition is estimated as (| sin(2(^i))^ ~ 0.23j)^. If we take a value of 
p=38 Debye corresponding to the full center to center distance the calculated intensity of 330De6ye^ 
is too large by a factor of about 3 to 4. The quantum calculations, however, reduce p for the cationic 
states by the rearrangement of the electron system to a value of about 20 Debye and the calculated 
intensity (dODebye^) is close to the experimental data of SODebye^ for viridis. From the experimental 

transition energies we find approximately 2w(if + f/^ ss 2700cm~^. Together with ^ f 1 , ^^ ^ 

sin^(0i) ~ 0.35 we estimate the coupling and energy difference as Ui = 1300cr7T,^^ and 2di = 800cm^^. 
The value of Ui is larger than the calculated 800cm^^. This discrepancy could be partly due to a 
geometry relaxation in the cationic state but also to the parameterization of the resonance integrals 
which was optimized with respect to transition energies. In our simulation (fig. la) we accounted 
for this shortcoming and used the larger value of 1300cm~^. If we assume that the intensity of the 
trip-doublets is mainly borrowed from the hole transfer transition the total intensity of the two trip- 
doublets is given by pUi2/{E2 — -Ei)^ ~ 25D^ for the calculated coupling of 1600c?Ti~^ and an energy 
difference of 7000c?Tt~^. It depends sensitively on the value of the resonance coupling Ui2- An increase 
of the coupling by a factor of 2 already gives more intensity for the trip-doublets than for the hole 
transfer band. The distribution of intensity over the two trip-doublet transitions depends on their 
degree of localization. In case of full localization of the lower trip-doublet on Pm the intensity ratio 
of 2:1 is essentially given by the localization of the lowest doublets. If on the other hand the two 
trip-doublets form a pair of completely delocalized states which might be a reasonable assumption as 
the resonance coupling U2 is larger in magnitude than Ui, the lower of the trip-doublet transitions 
gains about sin^(^i + f ) = 97% of the intensity and the higher transition becomes rather weak. The 
calculated splitting of the two trip-doublets is for the symmetric dimer 2U2 ~ 2000cm~^ but may be 
enlarged due to energy differences of the triplet states on the L and M halves. In our INDO calculations 
we found a rather small asymmetry only if we used the structure with the rotated acetyl on the M-half. 
Experimentally it is not possible to assign this state clearly which supports the conjecture of strong 
delocalization and the resulting low intensity. This conclusion differs from the assumptions made by 
Reimers and Hush [ 5 ] . If the trip-doublets borrow their intensity predominantly from the hole transfer 
band they will be essentially both polarized parallel to the vector p. This polarization is consistent 
with the experimental anisotropy of 0.298 it 0.028 [1]. From the quantum calculations we find a value 
of 0.27. We like to point to the fact that the magnitude and direction of p are largely influenced by 
the reorganisation of the electron system. For a transition parallel to the direction connecting the two 
magnesium atoms the calculated anisotropy has a much smaller value of 0.21. On the basis of this 
assumption Wynne et al [1] concluded that the intensity of the trip-doublet is not provided by the 
hole transfer. Instead they assume that the trip-doublet has the polarization of the Qy transition of 
Pm- The transition to the trip-doublet gains its intensity from the one-particle interaction between 
the two dimer halves as discussed by Reimers et al. A detailed analysis will be given in a forthcoming 
paper which shows that the transition dipole contains a contribution from the hole transfer transition 
and a smaller one from the Qy transitions of the two halves, the latter depending on the degree of 



localisation of the P^ states. 

Interpretation of the P* spectrum 

The electronic structure of the dimer in the Qy region is expected to show at least four states which are 
mixtures of two local excitations P£ and P^ and of two charge transfer states L^M~ and M~^L~ . The 
lowest excitation is mainly ascribed to the antisymmetric combination P£ — P^ with some admixture 
of the charge transfer states. This excitation carries most of the Qy oscillator strength as the transition 
dipoles of the two dimer halves are nearly antiparallel. The second excitonic component , the so called 
upper dimer band carries much less intensity and overlaps with the absorption bands of the accessory 
monomers. From the absence of a strong Stark effect in this region it was concluded [8] that the 
upper dimer band couples only weakly to the CT states. The CT states are predicted to be in the 
region between the upper dimer band and the Qx transitions of the dimer. They cannot be observed 
in the absorption spectrum as they have only little intensity and are inhomogeneously broadened by 
electrostatic interaction with the environment. In addition to these four basic Qy excitations we now 
incorporate an additional excitation in this frequency region which can be visualized as two triplet 
excitations on the two halves which are coupled to a singlet double excitation of the dimer. It might 
be identified with the weak band at 15800cm~^ denoted by a star in fig. 4 of [ 1 ]. According to our 
INDO results this double-triplet state P** couples strongly to the charge transfer states and may be 
important for the interpretation of the Stark effect of the dimer absorption. The coupling involves 
the transfer of an electron from an occupied orbital to an unoccupied orbital on the other dimer half 
and is therefore closely related to the coupling of the trip-doublets to the lowest cationic states. The 
energy of this state is underestimated by the INDO method for the applied parameters. As the lowest 
calculated triplet of the Bchl molecules comes out too low by about 1500cm~^ and the calculated 
energy of P** is close to P* or even slightly below [2,9,10 ] we expect the correct position of P** 
somewhat less than 3000cm~^ above P* which places it at the position of the lower and sharper P* 
absorption band of Wynne et al [1]. In the following we want to show that further support for such 
an assignment comes from the prediction of its intensity, its width and its Stark effect. We consider 
coupling of the following basic excitations: P£,P^, L^M^ , M^L^ ,P** described by the interaction 
matrix [ 2] 
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The factor w | in front of the HOMO-LUMO cross coupling matrix elements results from the proper 
spin multiplicities. Contrary to P^ the excited state P* is highly delocalized. So we start our analysis 
with a symmetric dimer with P(P£) = -E(P^) = E{*) and E{L+M-) = E{M+L-) = E{CT). The 
excitations can be classified as symmetric or antisymmetric with respect to the C2 symmetry operation. 
The interaction matrix reduces then to two blocks 
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with the new couphng matrixelements t/_ = (3l*m* + Plm ^+ = (3l*m* — Plm Ud = -\/6/3l*m 

and the zero order excitonic energies E{P*{-^)) = E{*) it W. 

(5) 
As the integrals (3lm and (3l*m* involve similar atomic overlaps in the ring-I region they have the 
same sign. Therefore the coupling between excitonic and CT excitation is enhanced for the negative 
states whereas it is reduced for the positive ones , i.e. the upper dimer band will be largely decoupled 
from the two other positive excitations. In addition our INDO calculations gave a strong dynamic 
correlation effect for the upper dimer band largely reducing the coupling to the CT states. The 
coupling matrix elements of the other states U^ and Ud show only minor changes as can be seen from 
table 2. 

Since the coupling [/+ is so small we treat it as a perturbation. Diagonalizing the strongly coupled 
states we get as approximate eigenstates of the symmetric dimer: 
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Substituting the experimentally identified energies and the calculated couplings gives the ordering 
s^ < s^ < c^ < c^ . The intensities for transitions from the lowest excited state P*(— ) are in this 
approximation 

/,(P*(_) ^ p*(+)) = - (Api - Apm) 

+ \e{P*(+))^E{CT(+)) + £;(p*(+))-£(P")j S-Pct 

/,(P*(_) ^ c'T(+)) = c+s_p,t - ^^p,,^^^±^^^^^^^\{^pL - Apm) 

/i(P*(-) ^ CT{-)) = ic_s_(ApL + Apm) 

As the dipole changes Ap of the isolated chromophores are small, the main contribution of the first 
three transitions comes from the large dipole moment of the charge transfer states pcT which from 
the INDO results has a magnitude of 32 Debye. So we predict for a nearly symmetric dimer that only 
the transitions to CT(+) and P** should carry significant intensities. 

Symmetry breaking interactions and Stark effect 

Breakage of the C2 symmetry causes a splitting of the two charge transfer states L^M^ and M^L^ . 
Such perturbations easily arise from different orientation of the acetyl groups or from asymmetric 
interactions with the surrounding , including the local electric field. We model the splitting of the two 
local CT states by a perturbation operator of the form 

H' = |(|L+M- ><L+M-| - |M+L- ><M+L-|) 

= |(c+s_|P*(-)><CT(+)| + c+c_|Cr(-)><CT(+)| (9) 

+s_s+|P*(-) >< P**\ + c_s+|Cr(-) >< P**\ + h.c.) + ... 

Here 5 is the perturbation parameter which admixes the states of different symmetry. The asym- 
metry does not effect in lowest order the intensities of the first three transitions in eq. 8. Only the 
transition to the fourth state CT(-) will borrow intensity from CT(+) as the two charge transfer states 
decouple into the localized states L^M~ and M^L~ . On the basis of these results we interpret the 
P* spectrum in the following way: The transition to the upper band at ~ 1500cm~^ is difficult to 
assign since it is weak and overlaps with infrared active vibrations. The observed band at 2700cm~^ 
represents the transition to the double-triplet state which is coupled strongly to the CT states. The 
broader band at 5000cm~^ contains transitions to the two CT states which are split and broadened 
by their interaction with the environment. This interpretation is also consistent with the different 



widths of the two bands. Only a small part of the width of the CT states is transferred. Figure lb 
shows a simulation of the P* spectrum using the calculated couplings from table 2. The energies of 
the transitions had to be adjusted slightly to fit the maxima of the experimental spectrum. Using 
the calculated couplings we have s_ = 0.37 ands+ = 0.45 . This gives a charge transfer contribution 
of 14% to the P*{-) state and a total intensity of the P* spectrum of s'^p'^ ~ lAODebye^ which is 
already close to the experimental value of 200Debye^ . The simulation with the parameters from table 
3 agrees quite reasonably with the observed intensities (fig.l). Further support for this assignment 
comes from the Stark spectrum. 

The energy of the lower dimer band depends on the asymmetry 6. Using perturbation theory of 
fourth order with respect to 6 and second order in s_ and neglecting the zero order energy difference 
of the CT states we find 

E{S)=E{P*{-)) 

2) ^- [e{CT)-E(P'-{~)) +■5+ (s(P")-E(P*(-)) ~ i?(CT)-£(P*(-)))) ^ ' 
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Even at zero electrical field asymmetric nuclear motions and static structural asymmetries con- 
tribute to the splitting of the CT states. Therefore we write 6 as the sum of a field dependent and a 
constant term 

6 = So + 2p,tF (11) 

Expanding the energy as a function of the electric field we then have 

E{F, So) = E{0, 60) + p{5o)F + 0(^0)^2 + . . . 
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which shows the effective permanent dipole moment and polarizability of the lower dimer band 
For a symmetric dimer the dipole moment is zero and the polarizability is approximately given by 

2 Pet /1Q\ 

"0 - -'~E{CT)-E{P*{-)) ^^^) 

The dipole moment as well as the polarizability grow with increasing asymmetry. The higher 
order terms are small as long as the asymmetry 6 is small compared to the energy difference E{P**) — 
E{P*{—)). For the special pair we estimate E{P**) — E{P*{—)) = 3000cm~^. The asymmetry which is 
necessary to explain the experimental Stark effect [9] is of the same order. Hence the higher order terms 
are in fact important. The Stark effect is strongly enhanced via the interaction with the state p** . For 
maximum mixture of the P** and CT(+) states we have s\ = 1/2. For E{CT) — E{P*) = 6000cm^^ 



and E{P**) — E{P*) = 3000cw, ^ this leads to an increase of 50% both for the dipole change and 
the polar izabihty. For s^ = 0.14 the calculated dipole change is without this enhancement effect 
(0.73(5/1000cm-i) + {)M{5 /imQcm-^f)Debye which yields a value of b.bDehye for 5 = 4000cm~i. 
In our INDO calculations such a large splitting resulted only for the acetyl rotated structure. 
Interaction with the P** state increases the permanent dipole moment by 17% to GADebye. This is 
roughly the correct magnitude consistent with our earlier assignment. The analysis of the experimental 
spectra, however, is complicated by its shape which contains contributions from the first and second 
derivative. Figure Ic shows a simulation using the calculated couplings from table 2 and the same 
adjusted energies as in figure lb. The calculated Stark effect contains a significant contribution from 
the large polarizability of the dimer which shows up as a first derivative. We like to mention that 
our model also allows to calculate the higher order polarizabilities which have been investigated in 
recent experiments [ 11]. For the absorption changes proportional to the fourth power of the field 
our calculations predict a mixture of up to the fourth derivative of the absorption spectrum with the 
major contribution from the third derivative. These seems to be in qualitative agreement with the 
experiment. 

Conclusion 

Our model calculations strongly suggest that a double-triplet excitation P** of the dimer is in the 
energy region of the internal charge transfer states and is strongly coupled to them. We identify 
this state with the narrow band of the experimental P* spectrum at 2700cm~^. This interpretation 
explains the intensity and the relative small band width as compared to the broad band at higher 
energies which is typical for CT excitations. Furthermore mixture of P** and the CT states enhances 
the Stark effect which helps to explain the experimental data. The relatively strong coupling involves 
the transition of an electron from an occupied orbital located on one half to an unoccupied orbital on 
the other half. This coupling is largely analogous to the coupling of the lowest cationic state of -P"*" 
to a trip-doublet state which has been invoked to explain the experimental -P"*" spectrum. Our INDO 
calculations are based on a structure model where the acetyl of Pm is rotated by 180° . This way the 
proper splitting of the CT states could be obtained. 

Acknowledgement 

This work has been supported by the Deutsche Forschungsgemeinschaft (SFB 143 and SFB 533) . 

References 

1. K.Wynne, G.Haran, G.D.Reid, C.C.Moser, P.L.Dutton, R.M.Hochstrasser, J.Phys.Chem. 100(1996) 
5140-5148 

2. P.O.J.Scherer, S.F.Fischer, in: The Photosynthetic Bacterial Reaction Center Center II, eds. 
J.Breton and A.Vermglio, Plenum Press New York 1992, p 193-207 

3. J.Breton, E.Nabedryk, W.W.Parson , Biochemistry 31 (1992) 7503-7510 

4. W.W.Parson, E.Nabedryk, J.Breton, in:The Photosynthetic Bacterial Reaction Center II, eds. 
J.Breton and A.Vermglio, Plenum Press New York 1992, p. 79-88 

5. J.R.Reimers, N.S.Hush, J.Am.Chem.Soc. 117 (1995) 1302-1308 

6. M.Plato, W.Lubitz, F.Lendzian, K.Mbius, Isr.J.Chem. 28 (1988) 109 



7. J.M.Ohlson, J.Tmnk, J.L.Sutherland,, Biochemistry 24 (1985) 4495 

8. P.O.J.Scherer, S.F.Fischer, in : perspectives in photosynthesis, ed. J.Jortner and B. Puhnian, 
Kluwer academic pubhshers 1990, 361-370 

9. P.O.J.Scherer, S.F.Fischer, in: The Reaction Center of Photosynthetic Bacteria III, ed. M.E.Michel- 
Beyerle, Springer New York 1995, p. 89-104 

10. P.O.J.Scherer, C.Scharnagl, S.F.Fischer, Chem.Phys.Lett. 197 (1995)333-341 

11. K.Lao, L.J.Moore, H.Zhou, S.G.Boxer, J.Phys.Chem. 99 (1995) 496-500 

12. B.S.Hudson, B.E.Kohler, Chem.Phys.Lett. 14 (1972) 299 

Figure caption 

figure 1: simulated spectra from the exciton model using the parameters from table 3 

(a) the experimental P"*" spectrum from [3] is compared with the calculated transitions shown as 
bars. As compared to the INDO results the HOMO-HOMO coupling has been enlarged by a factor of 
1.8 to reproduce the position of the hole transfer band and the doub-triplets are shifted by 2800cr?T,~^. 

(b) the calculated P* spectrum (full curve) is compared with a fit of the experimental data by 
two Gaussians (broken curve) as described in [1]. The scale is the same for both curves. The higher 
intensity seen experimentally for the lower band could be explained if the coupling of the CT states 
to P** and P*{-) were enlarged as compared to the calculated values. 

(c) the calculated Stark effect spectrum (curve) is shown together with the calculated absorption 
(reduced by a factor 0.05). The dimer band shows a mixture of first and second derivative with the 
minimum close to the position of the absorption maximum. The other excited states contribute only 
little to the Stark effect. The width of the dimer bands are as in earlier simulations of the room 
temperature spectra. The splitting of the CT states of 4000cm~^ is from INDO results on the M- 
acetyl rotated structure. The coupling of P** and the CT states is essentially the same as the coupling 
between the trip-doublets and the lowest cationic states of the dimer. The differences in table 3 are 
due to different reorganization effects. 



zero order 



Ui = -(3lm = 320cm-i 
U2 = \(iL*M* = -970cm"^ 
h2 = -\/^/3l*m = 565cm" ^ 
t/21 = -a/I 



U^2= y 
21 = -\/^(3lm* = W50cm-^ 



with CI 



Ui = 800cm- 1 
U2 = - 970cm" 1 
U12 = 1370cm- 1 

U21 = nrocm'^ 



Table 1: calculated coupling matrix elements of the lowest P^ excitations. The left column shows the CI 
matrixelements of the main configurations as a function of one electron resonance integrals (3 between 
localized orhitals of the neutral dimer groundstate. L(M) and L*(M*) denote the highest occupied and 
the lowest unoccupied orhitals of Pl{m) ■ The right column shows the effect of the reorganisation of the 
electron system in the cationic state. 





zero order 


with CI 


U- 


= Plm + Pl*m* = -2180cm"i 


U- = -2100cm"i 


u+ 


= Pl*m* - Plm = -1530cm"i 


U+ = - 202cm- ^ 


Ud = 


-\J\{I3lm* + f3L*M = 1690cm"i 


Ud = 1370cm- 1 



Table 2: calculated coupling matrix elements for the lowest excited singlet states of the dimer. The 
left column shows the one electron resonance integrals of the main configurations. The right column 
shows the dynamic correlation effect. 



Table 3: parameters of the simulation fig. 1 
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